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ABSTRACT 


The ability to predict surface currents can have a beneficial impact in several 
activities, such as Search and Rescue and Oil Spill Response, as well as others more 
purely scientific, operational or economic. The Naval Postgraduate School, in 
conjunction with the Romberg Tiburon Center and the University of California Santa 
Barbara, has been studying this purpose for the Coastal Response Research Center, 
University of New Hampshire and NOAA. So far the prediction was based on tide and 
persistence of the reminiscent current. Faced with increasing error under changing 
environmental conditions, further study of other influences became fundamental, in order 
to increase reliability. This study is a part of that effort by studying the impact of the 
wind-induced currents on forecasting. 

Based on a year and a half of wind and HF surface current readings, the wind 
surface current interaction is analyzed and quantified. Then that influence is plugged into 
the forecast algorithm. The final results show that the wind driven surface current is 
about 2% of the wind magnitude rotating around 50° clockwise, with coherence after 17h. 
The wind introduction into the forecast improved the accuracy, but only by an average of 
10%. The error still climbs with the variability of the environment, but knowing the wind 
influence allows other factors’ influences to be observed more accurately, such as the 
magnitude of the current itself. Forecasting is now done with 0.15 m/s plus or minus 0.1 
m/s at 95% confidence. 
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I. THE SURFACE CURRENTS FORECAST 

A. APPLICATIONS 

The use of the sea is not natural to man. Instead of moving with a solid surface 
underneath his feet and an almost intangible fluid surrounding him, at sea he has to deal 
with a much more viscous fluid. One of the major differences that this produces is that on 
earth, if he stops moving his body, or his car, he stops. At sea, he doe not. In order to take 
the best advantage of this distinction, understanding how this movement works is 
necessary. In the fifteenth and sixteenth centuries, the Portuguese voyages of discovery 
that navigated the entire globe in simple sailing ships relayed much of this knowledge. 
Charts of currents were drawn and pilots were written in each trip. 

With the passage of four centuries, several approaches have evolved to face the 
challenges of today. One of them is the remote evaluation of the surface current by HF 
RADAR (explained more in detail on II.B.2.a.). This method is excellent in filling the 
gap between satellite remote sensing measurements and the shore; however, it is not good 
for the Surf zone, for which other methods can be applied, like the one described in Puleo 
et ah, 2002. Depending on the frequency, HF radar can be used to read from around 3 to 
200 km off the coast. Surface currents are obtained with regular readings through this 
means, allowing the creation of a real time surface picture. Several works are still 
debating the accuracy of this system and what influences it. This study, focused on the 
forecast of the next 24 hours (24h) surface currents, takes the next logical step. Extracting 
currents using electromagnetic radiation (EMR) on the HE band implies that all the 
readings made are of the first two meters (2m) of depth, because EMR does not reach any 
further down and in this shallow layer, wind is considered to be a major driver of the 
currents. 

There are not many studies considering the interaction between the previously 
mentioned surface layer currents and wind, being almost all based on short time series 
and with limited resolution, most under the inertial period (24h). The existence of a 
database of two time series of both wind and current, with lengths over one year and with 
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an hourly resolution, presents an opportunity to unveil some properties of that interaetion 
that are currently poorly known. Based on these time series data (described in II.B.), this 
thesis’s aim is to first find the surface current component that is driven by the wind, by 
constructing a transfer function between them (described in Section II). This transfer 
function will then be used as an extra input for an already existing forecast model that 
relies solely on persistence and tides (described in I.B.), expecting to bring better 
accuracy and reliability to it. The results obtained for the wind/surface current interaction 
in this long time series may also be applied to many other future studies. 

Due to its practical implications, several institutions are interested in systems that 
allow this knowledge of the surface currents and in its forecast, because it has direct 
impact to every sea user that needs to plan any type of short term deployment. That is the 
case of the National Oceanic and Atmospheric Administration (NOAA) Coastal 
Response Research Center, who funded a study on “Delivery and Quality Assurance of 
Short-Term Trajectory Forecasts from HF Radar Observations,” of which this thesis is 
also a part. Better knowledge of the next 24h surface currents may be used to get better 
results in: 

1. Navy Type Operations 

Surface current forecasting has a significant impact on non-motorized activities, 
because propulsion reduces the dependency on the flow to which the user is exposed. 

a. Special Operations 

Deployment of divers or any type of raft maneuvering at the coast may 
benefit from a shortened mission time and energy, reducing risks and costs. 

b. Mine Warfare 

Locating the positions of drifting mines remains a modern reality because 
mines are still the cheapest and most effective way of making constraints to the use of a 
certain sea area or a harbor. To avoid or to sweep drifting mines, knowing where a field 
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is today will not help much tomorrow, because the mines drift with the current at their 
depth. If you are deploying, you need to know if the mine behavior will be what you 
intend. 


c. Flag Presence 

Most flag presence missions consist of long periods where a ship is 
constrained to a certain position or box. Taking advantage of the currents will result in 
better efficiency in the use of time and fuel. 

d. Maneuvers and Evolutions 

To maneuver a ship a long distance, the vector of her self propelled 
velocity must be added to the surface current velocity, also resulting in better efficiency 
in the use of time and fuel. 

e. Weapons Practice 

Target drift in time. 

In live surface fire practice, usually a target is deployed or towed at a very 
slow speed. More accurate knowledge of the current will predict target drift over time, 
allowing a better residence time in the exercise area with fewer maneuvers, making it 
safer. 


2. Coast Guard Type Operations 
a. Search and Rescue 

Most distressed platforms do not have, or only have partial, locomotion 
capability. Examples are man overboard, life saving rafts, debris or a damaged ship. 
Having the correct surface current scenario will allow two major advantages: 

(1) Reduces Position Uncertainty to the Distressed Position. 
The future target position is a function of time multiplied by the elements’ influence. 
When the rescue craft reaches the search zone, it usually has to face an area which size 
depends on the uncertainty of the surface current and the wind influence on the target 
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(and the initial position error). Redueing surfaee eurrent uneertainty will direetly reduee 
the seareh area, redueing the time needed to find the target, allowing better survivability 
odds. 

(2) Reseue Missions Savings in Hours and Fuel. Added to the 
advantage in (1) of faeing a shorter area, reseuers are able to plan a more aeeurate 
intercept course, arriving earlier than possible previously which results in an even tighter 
search area. On the scene, better sweeping plans can be made. 

b. Oil Spill Response 

Knowledge of the impact area, where a spill occurred and where it is 
going, allows the efficient deployment of the available means, reducing damage. For a 
good result, it must be coupled with the pollutant’s dispersive characteristics. Also, it 
helps track pollutants that may have come from tank washings. 

c. Illegal Traffic and Immigration 

Much of this traffic is done by small slow boats, for which the surface 
current is a determinant of the true course. Also, sometimes packages are sent adrift on 
purpose or to avoid capture. Knowing the currents will facilitate the collection of them. 

d. Patrol 

Better efficiency in the use of time, fuel and area covered. 

3. Fisheries 

a. For Those Who Monitor 

One typical periodic problem is the understanding of how a bloom will 
spread out. Knowing more about buoyant behaviors helps avoid the capture and 
consumption of organisms that may have toxic elements. Also, surface current 
knowledge gives a better understanding of the probable paths of buoyant (or slow 
velocity) organisms, such as eggs, larvae or youngsters. 
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b. 


For the Fisherman 


Coastal fishermen depend a lot on timing and their profits depend a lot on 
fuel eonsumption. The knowledge of how the eurrents will behave provides a tool to 
make transits and staying with fisheries’ spots more cost efficient. They also benefit from 
more accurate knowledge of the species behavior described in the previous paragraph. 

4. Meteorological/Oceanographic 

Other than the direct relevance, it is an important tool to models, as they need 
periodic assimilation of reality data to be accurate. This means can also control model 
forecasts’ quality by direct comparison with its own. 

5. Commercial and Recreational Navy 

Coastal transits can be planned to be more efficient in terms of fuel, time and 

safety. 


6. Habitat Impact 

Knowledge of buoyant behaviors makes it possible to understand the behavior of 
an atypical discharge from a river or sewer system, providing a tool to better manage the 
situation. 

B. ACTUAL MODELS 

Surface currents are usually collected by direct measurements using current 
meters or current profilers, or by remote sensing systems, such as coastal radars or 
satellites. In order to make forecasts, there is need of actual observations plus a model for 
future behavior, and forecasting ability will be in direct relation to the type of 
observations available. 

Retrieved data from satellite remote sensing is usually based on readings of sea 
surface temperature (SST), or other pixel differentiation information, to allow tracking of 
characteristics by maximum cross correlation. This gives resolutions good enough for 
coastal applications that range from 1 km to 25 km according to the frequency used, the 

latter for microwave, which is the only method able to penetrate clouds (Emery et al., 
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2005). They are then coupled with satellite altimetry information to create the 
geostrophic flow and make corrections, if a satellite is already on task. However, there 
are two primary disadvantages to using satellites: they are more expensive to operate and 
have less time resolution, despite the possibility of coupling with models in between. 
Since the time resolution is twice a day, a function of the orbit, it does not allow 
handling of smaller time scale phenomena. Those phenomena may not be very important 
in the open sea, but they gain substantial, if not major, influence near shore. One of such 
cases is the tides. Also, tides cannot be forecasted with poorer resolutions. No 
quantitative evaluation of the accuracy of these models was found in the references. 

The use of radar systems for this purpose is becoming widely used throughout the 
world because of the many unique applications that exist for real-time surface current 
mapping data. Also, it is a cheaper and less dependent system compared with the use of 
satellites, particularly if the user is interested primarily in coastal currents. Since the real 
time data is already collected, several academic, operational and commercial institutions 
have been trying to go further and forecast the surface currents. Following are 
descriptions of two models referencing this work. 

1. The Naval Postgraduate School Model 

This approach relies on the Coastal Ocean Dynamics Applications Radar 
(CODAR - a HF RADAR system) readings of surface currents, as described in II. B. 2. 
a.. These readings are then treated for Geometric Dilution of Precision (GDOP), which 
takes into account the uncertainty provoked by the angle of the azimuths between stations 
and reading position (the closer to 90 degrees, the better). The currents associated with 
tides are then removed by software from MathWorks called “t tide.” More information 
about “t tide” can be found in Pawlowicz et ah, 2002, which performs a best fitting 
frequency relation between the time series and the tidal harmonic constituents, taking the 
period of interest into account. This model uses a Cook, 2008, alteration to the software. 
After removing the tidal effects from the current, the remaining current is then averaged 
to the mean current for that day. Based on persistency, this current is then projected to the 
next 24 hours and is added to the tidal current that is again forecasted by the “t tide” 
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prediction software. The accuracy of this model is typically 7km distance after 24h, but 
frequently exceeds 20km.( Garfield et ah, 2007). This thesis builds on this model 
improvement by introduction of the wind influence. 

2. The Coast Guard / University of Connecticut Model 

Working with the U.S. Coast Guard, the University of Connecticut developed the 
Short Term Predictive System (STPS, Ullman et ah, 2003; O’Donnell et ah, 2005) for the 
U.S. Coast Guard Research & Development Center. This model also relies on CODAR 
readings treated for GDOP. In their work several essays were attempted. In the last, the 
surface current was divided into three parts: tidal flow, wind driven and mean. The tides 
and the wind influence are estimated by a different algorithm, the Gauss-Markov 
estimator, which divides what part of a signal is tidal or wind motivated and what is not, 
by comparison of time series, creating weights that are calculated from the 
autocovariance functions. This signal decomposition then allows prediction of the next 
25h by persistence, tide and predicted winds influence. One NDBC buoy is used for the 
wind data, with results more accurate than those from the previous Coast Guard results 
using only NOAA’s tidal data. Root Mean Square (RMS) errors were in the order of 10 
cm/s (ten centimeters per second) for each component of the velocity, u and v 
(corresponding to 14.14 cm/s for the vector). Forecasts for search boxes were then made 
with 20 cm/s RMS in order to cover uncertainty in all measurements and then compared 
with real data from drifters and statistical approaches by a Monte Carlo simulation. 

The main differences between the Coast Guard model and the one used in this 
thesis are: the length of the evaluation data set - seven days versus 1,5 years, the way the 
estimation for the tides and winds influences are calculated - Gauss Markov during one 
Month for the tides and during 1 day for the wind versus t tide for one Month and time 
series analysis for 1,5 years respectively, the wind data position (number of buoys) - one 
versus five different positions, and the different number of CODAR data stations - two 
versus 10 different stations. The STPS concludes that “...another analysis approach to 
quantifying the wind-driven circulation is called for.” 
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C. FORECAST AND THE INTRODUCTION OF WIND INFLUENCE 

All the works concerning forecast of surface currents point to the necessity of 
accurately knowing all its dependencies and their magnitudes. We know that the end 
result will be in direct proportion to the accuracy of the data and of the relations we plug 
in. One of those dependencies is the wind, as Ekman’s work proved by taking one step 
further our intuitive experience that the wind influences the movement of the water 
underneath it (Ekman, 1905). 

The work reported here was done after the initial results of the use of the Naval 
Postgraduate School model for the NOAA’s Coastal Response Research Center (Garfield 
et al, 2007), where the need to understand the wind’s influence on the surface current was 
identified as possibly relevant, among other factors. Despite this thesis’s focus on the 
wind, other dependencies will also be analyzed when doing error assessment. In order to 
reach a forecast model including the wind, a full time-series analysis study of wind / 
surface current interaction is made and shown in section II. After detailing this 
interaction, and having a transfer function created, a forecast model is put to work 
producing results, which will be described in Section III. Conclusions from these results 
are then extrapolated and discussed in Section IV. 
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II. THE WIND / SURFACE CURRENT INTERACTION 


A. PREVIOUS WORKS USED IN THIS STUDY 

Despite the long list of studies on air / sea interaction, few studies were found 
concerning strictly the interaction between the wind and the surface current, with 
quantified results for periods less than 24 hours. Nevertheless, portions of the following 
works contribute to the understanding of the wind surface current relation process: 

I. Ekman, 1905 

Although this thesis is based on mainly empirical results, Ekman provided the 
basic theoretical support for it (Ekman, 1905). Starting with the geostrophic equations 
with vertical friction and dividing the flow into geostrophic and frictional components, he 
related the wind stress and its influence on the oceanic flows as a function of depth, 
creating, among other results, what we now know as “Ekman Spiral” (Eigure 1). The 
Ekman theory assumes infinite ocean and constant eddy viscosity (Ay), which is an 
approximation of an already poorly known value (Radko, NFS class notes, 2007). 



Eigure 1. Ekman Spiral. - The angle between the wind and the surface current is 45“ 
clockwise in the Northern Hemisphere (Erom: http://oceanservice.noaa.gov, 
September 2008). 
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According to his work, the wind stress is given by: 


. dv 
oz 

where r is the wind stress and represents the water density as a eonstant. 


( 1 ) 


The amplitude in depth is 


f being the Coriolis term and z the depth. 
The direetion is shown by: 


( 2 ) 
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being the angle at the surface. 


For this study, the two main results are: 


• The angle of the eurrent at the surfaee, whieh is 45° to the right of the 
wind stress. 

• The amplitude of the veloeity shown by 

, -■ , \Ta' 


\ ^0£ I 
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(4) 


2. McNally, Luther and White, 1988 

This is a study on subinertial frequeney response of wind driven eurrents in the 
mixed layer. The data set eonsisted of 63 traeks of drifters, between 1976 and 1981, with 
results being binned at 5 days periods. The drifters operated at a 30m depth. The study 
deseribed a proeess of the interaetion between wind and eurrents, despite the different 
objeetive, analyzed depth and periods. Among the results are deseriptions of the 
following effeets indireetly related to this work. 
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a. The angle 

The angle was found to be around 25° to the right of the wind stress (at 
30m), which means the existence of the clockwise rotation in depth, despite the angle, 
was smaller than the expected by Ekman theory. 

b. Correlation 

The correlation was found to be significant between near inertial to 16 
days. No relation was found under the inertial period. 

c. Major Axis 

This was found to be from 75° to the right of the wind stress at near 
inertial to 15 ° at the lower frequencies. 

d. Lags 

Lags were observed from 30 ° at near inertial to zero, corresponding to 
under four hours in time. 

3. Foster, 1993 

This thesis studied the relation between diurnal surface winds and currents in 
Monterey Bay. The author analyzed the daily wind surface relation for September 1992 
in Monterey Bay, using two CODAR sites for the currents, two moorings and three 
coastal stations for the winds. The result was data on the behavior of the wind and current 
circulation and interaction in this zone of the Monterey Bay. The major results that 
impact this study are the current response to the breezes be under two hours (the time 
resolution of the data) and the fact that the more offshore currents rotate clockwise in 
time, to the right of the wind. 
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4. Xu and Bowen, 1993; Monismith et al., 2006 

This is a work on wave and wind driven flow interaction. From this study, the 
major contribution is the description of an Eulerian return flow (Hasselman’s drift) that 
cancels the Stokes drift, allowing for a measure of the uncertainty related to the waves’ 
influence on the currents. 

5. Yelland and Taylor, 1995 

This is a work on wind stress measurements from the open ocean. The values of 
the drag coefficient {Cd) for the theoretical approach made in II.D. are used according to 
this study’s conclusions; 

• For winds speeds at 10 m (Uio) from 3 to 6 m/s 

lOOOC^ =0.6 + 0.07t/io (5) 

• For winds speeds at 10 m from 6 to 26 m/s 

lOOOC^ =0.6 + 0.07t/io (6) 

6. Ardhuin, Chapron and Elfouhaily, 2003 

This is a work about the interaction between wind and waves. The result of this 
study that most contributes to this thesis is that for growing waves in short fetches, there 
is a wave induced stress that opposes the wind stress, consuming up to 10% of the 
energy. This directly influences the transfer function and the forecast that, practically, is 
only based on actual currents and wind forecast. To remove this influence, a measure of 
the wave height difference due to wind presence should be introduced - which it self is a 
function of magnitude and direction of the wind and of the waves height, wavelength and 
direction and of the fetch. 

7. Ullman et al., 2003 

These authors prepared a for the United States Coast Guard, Department of 
Homeland Security, on the use of CODAR technology in Search and Rescue (SAR) 
planning. This work was continued by the next study referenced. 
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8. O’Donnall et al., 2005 

This is also a study prepared for the United States Coast Guard, Department of 
Homeland Seeurity, on the integration of CODAR and a Short-Term Predietive System 
(STPS) surfaee eurrent estimate into the Coast Guard’s SAR Optimal Planning System 
(SAROPS). This study is a eatalyst and at the same time a measure of eomparison for the 
results of this thesis. As deseribed in its goals and approaeh were the same, but 

used more limited data and a different mathematieal analysis. One of its eonelusions 
mentioned the necessity of a work such as this thesis. This thesis results ended up being 
similar. 


9. Garfield, Paduan and Ohlmann, 2007 

A report prepared for the Coastal Response Research Center, this is a study 
funded by NOAA and the University of New Hampshire. One of the focuses of this 
report was to test the accuracy of the Naval Postgraduate School model, and it pointed to 
the necessity of studying the other interactions, starting with the wind. For this purpose, a 
comparison is made between CODAR forecasted positions and CODAR reading 
positions of simulated drifters after 24h. Its conclusions mention that “one environmental 
factor that is expected to play a role is wind,” despite inconclusive results. The study also 
mentions that there is “some connection between strong currents (...) and large RMS 
separation.” 

This thesis is an active part of the continuation of this work. 

10. Ardhuin et al., October 2008 

This study came to the attention of the author during the writing stage of this 
thesis. Coincidentally, it concerned exactly the focus of section II; the inducing of 
currents by the wind. Its conclusions bolster the ones reached in this thesis, having points 
in common such as the magnitude of the wind driven current measured around 2% and 
the angle between the velocity vectors observed from 10 to 40 degrees, based on a two- 
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year time-series analysis. Moreover, he study goes further in theorizing about the wave’s 
influenee, one of this thesis’s proposals for necessary future work , as pointed to by 
II.A.4. and II.A.6. 

B. DATA 

As mentioned earlier, this work is one of the first able to study two long reliable 
time series, with continuous one hour measures for more than one year. All data analyzed 
was gathered from June 1, 2006 to December 31, 2007. Figure 2 shows the time intervals 
for each station. This allows a wide spectrum of possibilities and test scenarios in space 
and time, as well as providing outputs of the time series analysis with a very good 
resolution, maintaining a continuous information flow with a high confidence level. 
Current data is gathered from 10 CODAR stations, with calculations of currents for five 
NDBC buoys positions, from where the wind data are obtained (Figure 3). 

Time Series intervals of matched Wind and Current data 



Figure 2. Matching data of Wind and Current - Each jump in the line corresponds to 
an interval with no data. 
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Figures. Current and Wind reading stations - Dots represent one example hourly 
coverage of the Pescadero CODAR station. (From: ArcGIS. Coastline from 
NOAA, National Ocean Service, 1994.) 

1. Winds 

Winds were obtained from NOAA’s National Data Buoy Center (NDBC; Figure 
4, Table 1), and from the Monterey Bay Aquarium Research Institute (MBARI; Figure 5, 
Table 2) moored buoys. 
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a. 


NDBC Buoys 



Figure 4. NDBC type buoy. (From: National Data Buoy Center, 
http://www.ndbc.noaa.gov , October 2008). 


Data from the following NDBC buoys / stations is used: 


# 

Name 

Distance to 
shore 

Type 

Position 

Elevation 

Anemometer 

height 

Water 

depth 

Watch 
circle radius 

46012 

Half Moon Bay 
(HMB) 

4NM 

3 m 
discus 

37.361N 

122.881W 

sea 

level 

5 m 

213.1 m 

170 yds 

46026 

San Francisco 
(SFR) 

18 NM 

3 m 
discus 

37.759N 

122.833W 

sea 

level 

5 m 

52.1 m 

127 yds 

46042 

Monterey 

(MRY) 

27 NM 

3 m 
discus 

36.753N 

122.423W 

sea 

level 

5 m 

2115m 

2175 yds 


Table 1. NDBC Buoys (From: National Data Buoy Center, 

http://www.ndbc.noaa.gov , October 2008.) 
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b. MBARI Buoys 



Figure 5. MBARI MO and Ml buoys. (From: http://www.mbari.org. October 2008). 
Data from the following MBARI buoys / stations is gathered from: 


Name 

Position 

Elevation 

Anemometer 

height 

Water 

depth 

MO 

6.835N 

121.899W 

sea 

level 

4 m 

70 m 

Ml 

36.750N 

122.020W 

sea 

level 

4 m 

1000 m 


Table 2. MBARI Buoys (From: Monterey Bay Aquarium Research Institute, 
http://www.mbari.org. October 2008; Depths (From: ICON, 

http://www.oc.nps.edu/~icon/moorings. November 2008). 

Notice that since MO and Ml are well inside the Monterey Bay (Figure 5), 
their data differ from the other’s positions because they are under the control of the 
typical Bay circulation and winds, of high topographic and land influence. 

c. Resolution and Accuracy 

The resolution of the anemometers is unknown but the readings are valued 
to 0.0001 m/s. The accuracy is also unknown, but considering the readings’ detail, 
variance is considered to be small. 
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d. Handling 

In order to enable treatment and comparison among different equipment 
readings, all the data must be formatted into a similar shape. Therefore, all of the wind 
time series are; 

• Hourly averages of the available 10 minute readings, from the exact hour - 
30 min to + 30 min; 

• Arranged in the form of a three dimensional vector of time (in hours), u 
velocity (seasonal, in m/s) and v velocity (meridional, in m/s); 

• Interpolated for gaps under 6 hours; 

• arranged so that seasonality, monthly, synoptic and daily influences are 
studied by subdivisions of the series in intervals of three months, one 
month, six days and one day. 

• calibrated such that spatial influence is studied by the different locations. 

e. Other Series Created 

A different time-series is also created for the wind readings of the 
Monterey and Half Moon Bay NDBC buoys, where the y axis is rotated for the winds’ 
largest variance azimuth, maintaining the x axis in perpendicular alignment. This will 
serve to interpret the angle between the two major variances of wind and current, 
extracting conclusions on interaction both for rotated axis for wind biggest variation only 
and for wind data rotated for the wind biggest variation and current to the current biggest 
variation. 

2. Currents 

Currents were obtained from all the in-range CODAR stations, for all the 
positions of the wind stations mentioned in the previous section from the buoys positions. 
A photograph of a typical CODAR transmit and receive antenna configuration is shown 
in Figure 6. Despite the increasing number of CODAR stations around the world, there 
have not been many opportunities to operate for a long period with so many stations in a 
single study, increasing the number of radial readings for each position. A short 
explanation of how CODAR technology appeared and operates is given in paragraph a. 
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a. 


CODAR Technology 



Figure 6. Point Sur CODAR station - The receiver antenna is closer and the 
transmitter further down. (From; CODAR, http://www.codaros.com , 
October 2008). 

CODAR consists of radar in the FIF band, which is the band that 
corresponds to wavelengths of the same magnitude of the sea waves. The first 
correlations between radar and waves appeared accidentally at the end of the Second 
World War, when echoes appeared in near shore radar stations. In 1955, “Crombie 
discovers experimentally that Bragg scatter produces strong HF sea echo” (CODAR, 
http://www.codaros.com , October 2008). 


CODAR uses the Bragg Scattering phenomenon (Figures 7 and 8) to 
measure the wavelength of the sea waves that are approaching or going away, in a radial 
direction, to each station. Having the wavelength, we calculate the speed a wave should 
have from the deep water dispersion relation. The difference of where the received 
frequency peak is and where it should be in that particular wavelength is a measure of the 
radial surface current by (Doppler, 1803-1853): 


/’ = /o 


1 + 


(7) 
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/' being the new peak frequeney, the expeeted peak with no current, v 
the radial current velocity and Cg the theoretical velocity of a wave with that wavelength 

in deep water. With two of these measurements, if properly intercepting, we obtain the 
two components of the velocity vector. 



Figure 7. Bragg scattering phenomenon - The biggest return in scattered energy is 14 
X (wavelength) of the sea wave. (From; Martin, 2004 in “An introduction to 
Ocean Remote Sensing,” Cambridge.) 
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Figure 8. Bragg scattering phenomenon. Received power by frequency - The 
difference of where the peak is and where it should be for that particular X , 
is proportional to the radial current velocity. (From; Paduan & Graber, 
1997.) 

For the frequencies used in these stations, the measurements are related to 
the first one meter of the water column; horizontally, the result for a single observing bin 
is based on backscatter from, approximately, 3x3 km of surface area. 

h. CODAR Stations 

All stations but the one off Moss Landing consisted of a standard 11.5-14 
MHz apparatus; Moss Landing operated in the range 24-27 MHz. Their approximate 
positions (in Geographic Coordinate System) are given in Table 3. 
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CODAR station 

Position 

Commonwheel 

37.911833N 

122.728167W 

Fort Funston 

37.712050N 

122.5013W 

Montara 

37.533717N 

122.519217W 

Pescadero 

37.252483N 

122.416567W 

Santa Cruz 

36.949217N 

122.0661W 

Moss Landing 

36.803667N 

121.787917W 

Point Pinus 

36.636783N 

121.9536W 

Naval Postgraduate School 

36.603333N 

121.87195W 

Granite Canyon 

36.439450N 

121.922217W 

Point Sur 

36.305767N 

121.90125W 


Table 3. CODAR stations (From: http://www.cencalcurrents.com. September 2008.) 
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c. Resolution and Accuracy 

The resolution for these stations is 0.04 m/s (from 
http://www.codaros.com. aeoessed in Oetober 2008). The most reeent studies show aeouraey 
from 8 to 18 em/s (from Wright, 2008, NFS Masters Thesis). 

d. Handling 

In order to be able to work with these time series all the data have to be 
arranged in a similar fashion. All of the eurrent time series are: 

• Hourly averages from the exaet hour -30 min to + 30 min, interpolated 
from the nearest valued positions to the buoys position. This is the way the 
data is presented by the CODAR system. 

• Arranged in the form of a three dimensional veetor of time (in hours), u 
veloeity (seasonal, in m/s) and v velocity (meridional, in m/s). This 
arrangement allows easy handling. 

• Interpolated for gaps under 6 hours. This allows a good continuity and 
guarantees no hiding of important features in the changing of the current. 

e. Other Series Created 

Other time-series were also created for the current readings of the position 
of the Monterey and Half Moon Bay NDBC buoys: 

• Where the y axis is rotated for the wind’s largest variance azimuth, 
maintaining a perpendicular x axis; 

• Where the y axis is rotated for the current’s largest variance azimuth, 
maintaining a perpendicular x axis. (Comparison of the results from steps 
(1) and (2) also provides information about the most common angle 
between surface current and wind.) 

• All the series created for the Monterey buoy have two more associated 
scenarios, where the tide influence from the current is tentatively removed. 
For this purpose the software t_tide is used (Pawlowicz, Beardsley and 
Lentz, 2002 with a Cook, 2008, alteration to the software). In the first, all 
of the tide constituents available in the software are used, once this time- 
series is long enough (in this case, one and a half years). In the second 
one, all but the SI constituents are used, in order to avoid removing the 
sea breeze influence. In both we have to take into account that the energy 
removed may not be just tide related, as occurs in periods of the order of 
the rotation of the planet’s rotation frequency. 
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• Seasonality, monthly, synoptic and daily influences are studied by 
subdivisions of the series in intervals of three months, one month, six days 
and one day. 

• Spatial influence is studied by the different locations. 

3. Both Wind and Current 

The time series of both wind and current are then cut in order to create time 
intervals where they both exist, to make comparison possible. The resulting matched 
intervals are displayed in Figure 2. 

From now on in this thesis, for reference simplicity, when a station is mentioned 
by place name it refers to the buoy position; e.g., “Monterey Station” means the data 
obtained for the position of the 46042 NDBC buoy. Because the Monterey station has a 
broader continuous time series, it allows us to extrapolate more results with a better 
resolution. That particular location is where most studies and essays are made, the others 
serving as comparisons. 

C. EMPIRICAL APPROACH 

Two independent paths are taken. The first starts with real data and tries to 
extrapolate relations from it. This approach relies on pure Time Series Analysis. The 
second, described in II D., is based on theoretical knowledge and serves as comparison to 
the empirical results. 

All of the following methods are applied in all of the time series described in B. 

I. Principal Axis (Biggest Variance) 

The objective is to have a measure of the axis which experiences the largest 
variance and its angle of measurement, for both wind and current. The result is plotted as 
Figure 9 and as a table in Table 4, where the two major semi-axes in the along vector 
direction are shown in full (using Bahr, 2007, software): 
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Elipses axis from 01-Jul-2006 + 1.5 years ; alfa=-58.5547deg 



U 

Figure 9. Principal axis analysis - Values correspond to variance and alfa is the angle 
between the two full lines. (Monterey station.) 

The table shows values for each term; 

Principal axis for W/C: from Ol-Jul-2006 + 1.5 years 



theta 

Maj. 

min. 

wind 

-30.4301 

4.8623 

1.9311 

current 

28.1246 

0.13539 

0.13232 


Table 4. Principal axis analysis. - Values correspond to magnitudes and angle 
between the major axis and north. (Monterey station.) 


The results of all the time series show that magnitudes of the variance keep within 
the same order. The angles vary more, getting dispersed as we reduce the time series 
length. The average angle between the two major variance directions is around 50 
degrees in the longer time-series as shown in Table 5. 
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Examples of results from the different stations in different periods 


Station 

Period 

Angle between major axis (°) 

MRY 

01JUF06 .. 31DEC07 

-59 

HMB 

01JUE06 .. 23NOV06 

-14 


20FEB07 .. 12JUF07 

-8 


04SEP07 .. 31DEC07 

-10 

SFR 

29NOV06 .. 07JAN07 

-2 


01MAR07 .. 26MAY07 

-15 


06SEP07 .. 31DEC07 

-19 

MO 

110CT06..28MAR07 

-39 


29MAR07 .. 25JUN07 

+21 


25JUN07 .. 31DEC07 

+25 

Ml 

16NOV06 .. 05JAN07 

-43 


24APR07 .. 24JUN07 

+1 


25JUN07 .. 30AUG07 

-6 


Table 5. Prineipal axis analysis. - Spaee and time variability (wind - eurrent.) 

Despite the tendeney for the longer time series, the high variability does not allow 
a good estimation of a typieal angle between wind and eurrent, demonstrating that it may 
be more a funetion of time and plaee. 

2. Complex Correlation (Kundu, 1975) 

Complex eorrelation gives a measure of eorrelation between two eomplex time 
series. The magnitude of the result shifts between 0 and 1, 1 being if the signal is the 
same. The average phase angle gives an idea of the physieal angle between the two 
phenomena. Results for the Monterey station are plotted on Figure 10, and results for all 
stations are tabulated in Table 6. 
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Complex Correlation W/C MRY JUL06 to DEC07 


0.8 
0.6 
0.4 
0.2 

^ 0 
- 0.2 
- 0.4 
- 0.6 
- 0.8 
-1 

-1 - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 

Real 

Figure 10. Complex correlation - Modulus is a measure of correlation and alpha is the 
physical angle. (Monterey station.) 

The magnitudes vary from 0.2 to 0.7. The angles tend to remain on the second 
quadrant around 35 degrees, gaining variance with the reducing of the length of the time 
series. These results mean that they are correlated, and there is in average 35 degrees 
between the directions of the two. This is expected because the angle is supposed to 
increase in time, until Ekman’s 45 degrees. The results for MO buoy present the highest 
variability but also the smaller correlation, probably due to the variability of the winds 
and currents inside the Bay. An attempt was made to verify if the correlation would 
increase when lagging the Current time series, but it decreased, pointing to zero lag. 


Examples of results from the different stations in different periods 


Station 

Period 

Magnitude 

Angle n 

MRY 

01JUE06 .. 31DEC07 

0.49 

-32 

HMB 

01JUE06 .. 23NOV06 

0.31 

-56 


20EEB07 .. 12JUE07 

0.66 

-47 


04SEP07 .. 31DEC07 

0.48 

-56 

SER 

29NOV06 .. 07JAN07 

0.66 

-1 


01MAR07 .. 26MAY07 

0.81 

-16 


06SEP07 .. 31DEC07 

0.62 

-23 


T-1-1 I .1 I I 

- mod=0.48542 alfa=-31.7407 



J_^_ 1 _I_I_I_I_L 
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Station 

Period 

Magnitude 

Angle (“) 

MO 

110CT06..28MAR07 

0.34 

-26 


29MAR07 .. 25JUN07 

0.27 

-83 


25JUN07 .. 31DEC07 

0.17 

-45 

Ml 

16NOV06 .. 05JAN07 

0.5 

-34 


24APR07 .. 24JUN07 

0.67 

-17 


25JUN07 .. 30AUG07 

0.69 

-13 


Table 6. Complex eorrelation. - Spaee and time variability. 

3. Scatter Comparison 

Seatter plots give an idea of eorrelation and permit a mathematieal best-fit 
expression (e.g., Figure 11, Figure 12, Figure 13). A perfeet eorrelation happens when all 
the pairs of values (points) of the two time series design a mathematieal expression. This 
will allow elaboration of the best transfer function relation in terms of magnitude. The 
one presented here considers the full energy spectrum, and the best fit is accomplished 
with a quadratic expression. 

This analysis can also be done by frequency bands if the choice is to transfer by 
frequency. As we will see later, the degree of added complexity for such an approach will 
not increase significantly the results of the first one. 
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Figure 11. Scattering diagram for U and best fit expression. (Monterey station, Jul 06 
to Dec 07.) 


Scatter for W vs CMRY 
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Figure 12. Scattering diagram for V and best fit expression. (Monterey station, Jul 06 
to Dec 07.) 
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Scatter for W vs CMRY 
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Figure 13. Scattering diagram for Modulus and best fit expression. (Monterey station, 
Jul 06 to Dec 07.) 

After applying it to all the series, including the rotated ones, the best proxy for the 
typical expression relating near surface wind to surface current, in terms of magnitude of 
the vector, is: 

C = 0.018t/5+0.03 ( 8 ) 

where C is the surface current in m/s and U 5 is the wind at 5 m height, also in m/s. 

The quadratic term shown on the figures is always positive and negligible, with 
some exceptions in the modulus when it goes negative. The results presented here are 
restricted to the linear terms in the fit. Those values multiplying by U5 vary plus or minus 
0.02. But, overall, the results show surface current magnitudes are very close to 2% of the 
wind speed magnitudes (Table 7). The choice of 0.018 is supported by the error analysis 
and tuning done in the forecast section. The last constant varies plus or minus 0.05. This 
variation tends to be smaller in winter and fall and in the fits based on longer time series. 
This expression will be tuned in the predicting algorithm in Section 111. 
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Examples of results from the different stations in different periods 


Station 

Period 

U 

V 

Modulus 

MRY 

01JUE06 .. 31DEC07 

0.006X-0.028 

0.018X+0.018 

0.007X+0.128 

HMB 

01JUE06 .. 23NOV06 

O.OOlX-0.063 

0.021X+0.096 

-0.016X+0.217 


20EEB07 .. 12JUE07 

O.OOOX-0.047 

0.017X-0.068 

O.OOlX+0.178 


04SEP07 .. 31DEC07 

-0.006X-0.045 

0.018X-0.013 

0.007X+0.167 

SER 

29NOV06 .. 07JAN07 

0.013X-0.031 

0.022X+0.037 

-0.016X+0.162 


01MAR07 .. 26MAY07 

0.013X+0.000 

0.018X-0.033 

0.006X+0.175 


06SEP07 .. 31DEC07 

0.007X-0.012 

0.023X+0.004 

0.013X+0.106 

MO 

11OCT06 .. 28MAR07 

0.005X-0.014 

0.007X+0.001 

O.OOlX+0.099 


29MAR07 .. 25JUN07 

O.OOlX-0.089 

0.004X+0.032 

0.006X+0.114 


25JUN07 .. 31DEC07 

0.004X-0.039 

0.006X+0.017 

0.006X+0.104 

Ml 

16NOV06 .. 05JAN07 

0.008X-0.039 

0.016X+0.041 

0.007X+0.111 


24APR07 .. 24JUN07 

-0.009X+0.025 

0.008X-0.028 

-0.005X+0.149 


25JUN07 .. 30AUG07 

0.005X+0.001 

O.OlOX+0.007 

-0.015X+0.177 


Table 7. Best fit expression. - Spaee and time variability. 

4. Auto Correlation (McMahan and Wyland, NFS Lab Software, 2007) 

Autoeorrelation of a given time series allows for a better understanding of the 


time seales involved. Eaeh oseillation represents the time seale of an aetive agent. 
Examples are shown for wind and eurrent in Eigure 14 and Eigure 15, respeetively. 
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auto correlation Wind U from OIOul-2006 + 24h 



X10 


auto correlation Wind V from 01 vJul-2006 + 24h 
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auto correlation Wind U from 01-001-2006 + 24h 



auto correlation Wind V from 01 ^ul-2006 + 24li 



t(h) 


Figure 14. Autocorrelation for Wind - The bottom is a zoom of the first positive 2000h 
3 Months. Daily, synoptic and monthly oscillations are present. (Monterey 
station, Jul 06 to Dec 07.) 
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Concerning the wind, the influence of the sea breeze is easily seen on the U 
component in the small oscillations of 24 h. On both, we can also see an oscillation 
around the 6 days (synoptic scale). In U, monthly and a seasonal oscillations are also 
present. 
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auto correlation Current V from OIUul-2006 + 24h 




1 1 


-1 -0.5 


0 

T(h) 


0.5 


1.5 

xIO* 


auto correlation Current U from OIUul-2006 24h 
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Figure 15. Autocorrelation for Current - The bottom is a zoom of the first positive 
2000h ~ 3 Months. Daily, synoptic and monthly oscillations are present. 
(Monterey station, Jul 06 to Dec 07.) 
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For the current, daily, synoptic and monthly oscillations are noticed as well. The 
tidal influence can be easily seen from the difference between Figure 15 (simple current 
measurement) and Figure 16 (after tide removal). Daily, synoptic and monthly 
oscillations are present, but the daily one is much less energetic. All the daily oscillations 
tend to be bigger in the summer. There is no important change looking between stations, 
except in the magnitude of the oscillations. 


a: 


Tides out - auto correlation Current U from 01>Jul-2006 + 24h 



Tides out - auto correlation Current V from 01 ■Jut-2006 + 24h 
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Figure 16. Autocorrelation for Current with tides removed - The bottom is a zoom of 
the first positive 2000h ~ 3 Months. Daily, synoptic and monthly 
oscillations are present. (Monterey station, Jul 06 to Dec 07.) 
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5. Power Spectra and Straight Coherence 

The power speetra of wind and eurrent give the distribution of energy versus 
frequeney (Figures 17 and 18). Based on this, coherence gives the significance of the 
relation of the signals for each frequency and the phase difference between them (Figure 



Figure 17. Power spectra of Wind and Current - Dashed at 95% confidence. (Monterey 
station, Jul 06 to Dec 07.) 

The peaks for the daily and 12h periods are easily spotted in both wind and 
current. Winds are especially energetic in the sea breeze period while the currents show 
almost equivalence on both, maybe because of the previously mentioned wind and tide 
correlations. After the tide signal removal attempt, both current peaks are slightly 
smaller, representing less energy. When removing tides by best fit of frequencies, we are 
removing also part of the wind influence which shares those same frequencies - the sea 
breeze. 
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Tides out - Gxmean Gymean MRY 



Figure 18. Power spectra of Wind and Current with tides removed - Dashed at 95% 
confidence (Monterey station, Jul 06 to Dec 07.) 

Adding the variances both with and without tides, we notice the slight decrease in 
energy that is removed. The results are summarized in Table 8. 


Difference of with and without tides 


wind variance 

19.5694 

tides in - current variance 

0.016719 

tides out - current variance 

0.015593 


Table 8. Added variances ~ energy. (Monterey station.) 

These series are significantly coherent for periods longer than a day. During this 
interval, the phase augments from near zero to around 60°. In the shorter time series the 
coherence is consistent only after 40 h, and the phase tends not to exceed 50“. When done 
by components (u and v) it is coherent after approximately 12 h. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

f (CPH) 


Figure 19. Coherence and Phase between wind and current - Dashed at 95% 
confidence. (Monterey station, Jul 06 to Dec 07.) 


6. Rotary Spectra (Gonnela, 1972 and Mooers, 1972, Bahr, 2007 
Software) 


Rotary spectra subdivides a signal into two rotating components, one clockwise 
and the other counter-clockwise. This allows one to analyze the coupling of two signals 
in terms of their rotating behavior, giving independence from the axis choice. Results for 
wind and current for the Monterey station are shown in Figure 20 and Figure 21, 
respectively. 



Figure 20. Rotary components for wind - Notice the two peaks for 12 and 24h. 
(Monterey station, Jul 06 to Dec 07.) 
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Figure 21. Rotary components for current, tides in and out - Notice a small decrease in 
energy. (Monterey station, Jul 06 to Dec 07.) 

The two diurnal and semi-diurnal peaks are distinct in both wind and current, 
reflecting the importance of such frequencies. Notice that the clockwise component is 
bigger, the peaks being more similar to the 24h point in the wind and the 12h in the 
current. The current peaks are only a little smaller after the tides’ removal. 

The “rotary coefficient” is a measure of the purity of the signal in terms of 
rotation; a fully rotating signal is 1, a non-rotating is 0. The rotary coefficient for the 
wind and current records at the Monterey station are shown in Figure 22 and Figure 23, 
respectively. Distinctively, for periods under 12h, the current is more rotational than the 
wind, which only shows two small peaks in the frequencies of interest. 
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Tides in - Rotary Coefficient for W MRY. From 01-Jui-2006 + 1.5 years 



Figure 22. Rotary coefficient for Wind - A fully rotational movement would be valued 
1. (Monterey station, Jul 06 to Dec 07.) 


Tides in - Rotary Coefficient for C MRY. From 01-Jui-2006 + 1.5 years 



Tides out - Rotary Coefficient for C MRY. From 01-Jui-2006 + 1.5 years 



Figure 23. Rotary coefficient for Current, tides in and out. (Monterey station, Jul 06 to 
Dec 07.) 
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7. Vector Cross Spectral Analysis (Mooers, 1972 and McNally et al, 
1989, After lessen and Bahr, 2007 software) 

Vector cross spectral analysis is a tool to check the coherence between two rotary 
signals. Plugging in wind and current, the coherence results for the Monterey station are 
shown in Figure 24. 



Phase 



-0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 

f (CPH) 


Tides in - Coherence for MRY from 01-Jul-2006 + 1.5 years 



Figure 24. Complex Coherence and Phase - Lower image is a zoom of the coherent 
part. (Monterey station, Jul 06 to Dec 07.) 

The application of this method allows one to calculate, besides the coherence, the 
angular shift in time. Table 9 shows the average results for the coherent-positive part. 
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Complex Coherence results 


Time (h) 

Angle (“) 

20 

30-35 

30 

40-45 

60 

45-60 


Table 9. Wind and Current shifting of angles in time. (Averages for all the stations 
and periods.) 


After 17h almost all the series are coherent and have a much slower variation. For 
periods between 0 and 17h the results were inconclusive. This treatment is done also to 
the velocity square, to search for a difference between pure wind velocity and wind stress 
(function of the velocity square). No significant difference was found between the 
coherence results for wind or wind stress. For the de-tided series, the relation keeps its 
results except during the 24 and 12h periods, where very narrow positive and negative 
peaks appear (Figure 25), probably a result of energy removal in those frequencies. 


Tides out - Coherence for MRY from 01-Jul-2006 + 1.5 years 



_0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 

f (CPH) 
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Tides out but S1 - Coherence for MRY from 01-Jul-2006 + 1.5 years 



Phase 



_0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 

f (CPH) 


Figure 25. Complex coherenee and Phase, tide removed - Top, all eonstituents. 

Bottom, all constituents removed but SI. (Monterey station, Jul 06 to Dec 
07.) 

Comparing results for other stations (not shown), the MO buoy is only coherent 
after 24h at a 50“ phase. The same can be said for Ml in the summer. 

The empirical approach gave us estimation of the interaction between the wind 
and the surface current in its magnitudes and phase shifts. This allows building a transfer 
function that will later be used to forecast the current, having the forecasted winds in 
consideration of various factors. 

D. THEORETICAL APPROACH 

In this approach, theory is used to compare with the empirical results. 
Calculations are made using Ekman’s works theory, and a function is created 
considering; 


Pa^D^ 10 ’ 

And at the same time 


(9) 
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~ Pa^V 


( 10 ) 


So, adding the two and rearranging, 

du = ^^-dz-U,,\ ( 11 ) 

Po 

The magnitude of Uio is taken to be approximate of U 5 , for it is the top part of the 
exponential wind growth in height; the Yelland and Taylor, 1995 drag coeffieient is used 
as funetion of the wind speed. 

E. COMPARISON 

Having the empirieal and the theoretical result transfer functions, and using Uio = 
10 m/s we obtain: 

1. From the results of the study in equation ( 8 ), 

du=0.21 m/s. 

2. From the theoretical result of (11), 

a. Ekman depth = 50m and Av=0.01 m /s du=0.68 m/s; 

b. Ekman depth = 10m and Av=0.01 m^/s du=0.27 m/s; 

c. Ekman depth = 50m and Av=0.3 m^/s du=0.02 m/s; 

d. Ekman depth = 10m and Av=0.3 m /s du=0.01 m/s. 

This means that the theoretical result is highly dependent on a good estimate of 
eddy viscosity, and that using a typical value (Ay-O.Ol), Ekman depth must be much 
shallower. Using typical values, Ekman theory overestimates the change in the surface 
current velocity. 

F. CONCLUSIONS 

Having studied the wind / surface current interaction in different perspectives, 
times and spaces, the relationship found can be described as follows. 
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1 . 


Transfer Function 


a. Magnitude 

The surface current generated is around 2% of the wind speed. 

b. Angle 

The angle measurement increases in an unknown way until around 17h. 
From then on, the angles vary from 35“ to 60“, according to Table 9. 

2. Notes 

Despite these results, the scatter plots show a large dispersion, which means that 
there is a high variability. 

The unknown growth until 17h includes wind wave interaction. It includes also 
components of transfer from 0 to 35“ of some magnitude. 

In this study, rotating the axis to biggest variance or long-shore, cross-shore, does 
not bring any significant change, probably because of the meridional characteristic of the 
Californian coast. Nevertheless, when using both wind and current axis to their respective 
biggest variances, the component minor wind axis versus minor current axis does tend to 
zero. 

The tidal removal does not bring major changes to the transfer function as a bulk. 
If transferred by frequency, care has to be taken especially in the inertial periods and its 
subdivisions. 
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III. THE FORECAST 


A. FORECAST WITH CODAR 

Forecasting with a remote sensing means like CODAR has the enormous 
advantages of independence and timing—independence because it does not rely on any 
other system to operate, and timing because it is updated every hour. This is more than 
enough to get good estimates for the surface currents variability, making data available 
when needed. Nevertheless, relying only on previous surface currents readings from our 
own system has its dangers and limitations. As will be shown, coupling our readings with 
those of other systems yields more reliable and accurate results. 

1. Reliability of a Forecast 

CODAR has proved to be a very reliable system, so far. Nevertheless, its results 
may be misleading. Once it uses electromagnetic radiation to operate, it is subject to all 
its difficulties and laws, as described by the following examples. 

• If an object with electromagnetic potential and length comparable to the 
transmitted wavelength (approximately 10m) goes near the antennas, the 
reception pattern changes, leading to deviation in the result; 

• The presence of a moving target, such as a ship, may introduce an error 
due to strong backscatter; 

• A change in the background radio frequency properties will influence the 
signal also. 

These and other problems have been addressed several times, and are adequately 
explored in several works (some of which can be found on http://www.codaros.com) . The 
use of a second system allows us to detect possible malfunctions by comparison, and 
gives a different set of solutions as backup. This is always preferable. 

2. Accuracy of a Forecast 

The accuracy of the system is dependent on the errors of CODAR (8 to 18 cm/s, 
from Wright, 2008), plus the ones associated with the estimation of the next 24h behavior 
(the assumptions and the errors of the systems we are building our prognostic on), plus 
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the error of what we eannot or do not prediet. If we introduee to our algorithm a good 
deseription of the future eonditions and have good relation funetions, we will obtain a 
more realistie result. 

So, eoupling with other systems, as upload or download, is advantageous. We ean 
obtain data from tides, winds, drifters, eurrent meters, and models—and ean assimilate our 
results baek into the models. 

B. PREVIOUS ATTEMPTS 

The previous attempt model is deseribed in I B. The results have been generally 
aeeurate for persistent weather eonditions, but tend to get worse with ehanging weather 
eonditions. 

1. Problems 

The model was based only in tidal and persistenee analysis. Onee any of the 
eonditions—wind, waves, or any other parameter—ehanges, the persistenee does not 
follow. In weather, signilieant ehanges ean happen in any 24h period. 

2. Results 

As pointed in I.B., there is still laek of eonfidenee in this model any time there is a 
signilieant ehange in the weather (namely in eurrent veloeity and possibly wind). A 
summary of results for the previous tidal model applied to the region offshore San 
Franeiseo is shown in Figure 26. Here a pixel grid of pseudo drifters release points was 
used to ereate predieted trajeetories for 24h, whieh was repeated eaeh hour for the period 
August 2006 through July 2007. The plot shows the rms spatial differenee between the 
end points of the predieted trajeetories eompared with the observed trajeetories one day 
after. The negative eorrelation with number of drifters (lower panel) refleets the faet that 
higher speed eonditions result in more pseudo drifters, leaving the sampling grid within a 
24h period. The direet effeet of wind speed on the spatial aeeuraey of the predieted 
trajeetories is less elear. Nonetheless, wind (speeifieally wind ehanges) is hypothesized to 
explain some error in foreeast eurrents. 
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Figure 26. Results for RMS eompared with wind and number of pseudo drifters - RMS 
for 24h. The number of pseudo drifters present may be used as a proxy for 
current speed. (Image from Garfield et al., 2007.) 

C. THE INTRODUCTION OF WIND 

One attempt to reduce forecast error is the introduction of wind influence. A study 
of the wind surface / current interaction was made for this thesis and is described in 
Section II. The objective was to build a transfer function that allows us to plug in the 
wind driven current into the forecast. Several approaches were considered: transfer by 
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frequency, transfer by different division of spatial and time components, transfer as a 
vector, transfer with tide analysis, and transfer based on persistence. The approaches that 
seemed to be more efficient for practical purposes (accuracy versus simplicity) narrowed 
to only one wind/current transfer function and transfer as a vector. Transfer by frequency 
was not chosen because of the relatively small magnitude of the end results. Tide removal 
and simple persistence conditions would have yielded similar results. The models tested 
were based on the diagram in Figure 27, regardless whether or not tides were removed. 

The winds had different methodologies applied also in terms of what time slice 
constituted future for computing wind differences and whether or not to use wind 
difference compared with removing and then reintroducing the full wind-driven current. 
The results were practically the same. In this hindcast study, forecast wind is taken from 
the actual readings, simulating a zero error wind forecast. The methodology tested here 
can, however, be immediately transferred to winds forecasted by an atmospheric model 
such as the Navy's COAMPS model. 



Figure 27. Basic flow chart diagram of the forecast procedure - Several adaptations 
were attempted. 

1. With Tides Treatment 

This attempt corresponds exactly to the previous diagram. The best results were 
obtained using the following procedure. 

We are at t=0 and we want the next 24h surface currents forecast, so we 
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• Take the last 24h eurrent readings; 

• Remove the tides by “t_tide” using all the previous month’s currents to the 
best fit. Conducting the detiding with using one month of data or just three 
days of data both returned good results as shown in Table 10. 

Detidal associated errors 


RMS (m/s) 

3 Days 

1 Month 

U 

0.1087 

0.1073 

V 

0.1027 

0.1145 

modulus 

0.1515 

0.1593 


Table 10. Detidal errors with 3 days and 1 Month data - The calculation of the Root 
Mean Square error will be discussed in III. E. (Monterey Station.) 

• Forecast the tide for the next 24h using “t_predic” and the structure of 
constituents that resulted from the detiding; 

• Low pass fdter the residual current to 18h; 

• Average the residual current; 

• Take the last 24h wind readings; 

• For each hour in the future h, get the forecast wind from that hour until 
t=0 and the wind readings before that until h-24; 

• Low pass filter both to 18h; 

• Do the average for each; 

• Take the difference of the winds; 

• Use the transfer function to get the current associated with the wind 
difference; 

• Sum the residual current, the tide and the wind current. 

Figure 28 is an example of the obtained result, in a day where a major wind 
change took place and the wind influence is noted. The resolution is hourly. 
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(m/s) (m/s) 


01/Jan u rms=0.060447 



01/Jan V rms=0.047992 



Figure 28. One day forecast with tide treatment - The wind changing its influence. 
(Monterey Station.) 
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2 . 


With Persistence 


The procedure sequence is the same as the above but with the major difference of 
not removing the tidal influence. This method uses the last 24h current readings as the 
current for the next 24, summed with the wind difference driven current. Example results 
using the persistence model are given in Figure 29. 


01/Jan u rms=0.057852 



01/Jan V rms=0.06595 



Figure 29. One day forecast with persistence - The wind changing its influence. 
(Monterey station.) 
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D. TUNING 


The wind interval to aeeount for the difference between forecast and real is 
chosen after several different runs bases on shifting the future wind interval for several 
periods forward. The one with least error is the one designated, where the wind interval 
ends at tp (Figure 30). Several values around the transfer function were also attempted. 
The best results are obtained with a 50“ clockwise rotation, and with a magnitude 

Q,„.=0.018t/,+0.018. ( 12 ) 



wind to take the difference 



last 24h wind 

1 

t 

t ^ 

p 


to 


Figure 30. The difference between winds - Several positions of the future wind interval 
were attempted for a given time of prediction (tp). 

E. COMPARISON 

RMS is calculated on a daily basis with and without introducing the wind . RMS 
is based on the difference between our calculated surface current for each hour and the 
actual CODAR reading at that time. This means that we are considering CODAR as 
“true.” We must take into account the CODAR intrinsic variability when analyzing 
results. The end RMS result for each component u and v is: 


with the total 


RMS 


component 


{reading - forecastY 


N-l 


(13) 


= ^RMS^^+RMSY . (14) 

The results are given in Table 11, considering a normal distribution. 
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(m/s) 

average 

Variance 

Standard 

deviation 

95% conf 

(+/-) 

99% conf 

(+/-) 

> 

RMS 

0.1593 

0.0022 

0.0471 

0.0923 

0.1213 

o 

a 

a> 

T3 

With wind 






RMS 

0.1667 

0.0028 

0.0529 

0.1036 

0.1363 

H 

No wind 






dj 

o 

RMS 

0.1427 

0.0023 

0.0482 

0.0945 

0.1242 

C 

fD 

With wind 






u 

(D 

RMS 

0.1495 

0.0028 

0.0528 

0.1035 

0.1360 

CLh 

No wind 







Table 11. Root Mean Square eomparison - Average daily values. (Monterey station, 
Jul 06 to Deo 07.) 


The statistioal results for the 18-month study period show that the wind influence 
is one order of magnitude inferior to the total error. Use of persistence or tide treatment 
yields practically similar results, with simple persistence of today’s observed low pass- 
filtered currents giving slightly better results. 

Looking at other station results, station MO gets the least error results, followed 
by Ml, probably due to the repetitive nature of the flow in that position, or perhaps 
because it is more protected from unaccounted variations. Results distributed by station 
are shown in Figure 31. There is a visible relation between error and stations, showing 
that the error causes are probably associated with the position and its typical patterns of 
wind, currents, the radar station itself, and other features. This will be studied more 
closely in Section III.F. Interestingly, the distributed results show clear, positive 
correlation between forecast error and the mean station current speed or current 
variability, which is consistent with the previous findings of Garfield et al., 2007. 
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Variability at 95% (m/s) Current magnitude (m/s) 
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Figure 31. Relation of RMS error and the stations - Top is magnitude of the current vs. 

error. Bottom is variability vs. error. Based on the tide plus mean with 
transfer function model forecasts. 
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Seasonality does not seem to play a major role, but the small quantity of data does 
not allow a definitive conclusion. The four-season results for the Monterey station are 
shown in Figure 32. 
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Figure 32. Relation of RMS error and seasons - Top is magnitude of the current vs. 

error. Bottom is variability vs. error. (For Monterey station, based on the 
tide plus mean wind transfer function forecasts.) 

The long records from the Monterey station can be used to illustrate the temporal 
behavior of the forecast error and the influence of the wind-based corrections. The daily 
RMS current speed errors are plotted in Figure 33 with and without the wind-based 
corrections. Almost everywhere, the wind-treated forecasts, blue and black, are lower. 
Red and green tips are easily seen constantly at the top, despite the small difference. 
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RMS compare - Wind influence 
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Figure 33. Zoom of a plot with the different calculated RMS - Wind influence, “p” is 
persistence, “nw” is no wind. (Monterey station.) 

Zooming in on portions of the record (Figure 34) allows us to see the nature of 
some types of errors. As one example, sometimes, not often, the RMS associated with the 
tides treatment forecasts goes up with no apparent cause, probably the result of a bad 
frequency fitting. 
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Figure 34. Zoom of a plot with the different caleulated RMS - Possible tide problem, 
“p” is persistenee, “nw” is no wind. (Monterey station.) 

F. ERROR CAUSES 

If we take into account the average current velocities in this area, 0.17m/s, errors 
of 0.15m/s represent almost the same magnitude. Nevertheless, considering the main 
objectives of the application, having a measure of the uncertainty is a good result. 

If we have a distress in a position where the rescue mean takes T time to fulfill the 
transit, with a C uncertainty current, after the transit the search area will be at least a 
circle of radius TxC = M meters. Reducing the C will proportionally reduce M, which 
means that it will take less time to find the target. 

The smaller the error, the bigger the odds of saving a life or reducing an oil spill’s 
impact. With this in mind, continuation of this research should be a priority of the naval 
sciences. 
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1. Accuracy Expectancy 

In order for us to have an idea of what to expect in terms of accuracy, we need to 
layer known uncertainties related to the surface currents. 

a. CODAR (From Wright, 2008) 

In a good scenario, we will have approximately 0.08m/s accuracy. 

b. Winds 

The influence is unknown, but considered to be small, according to our 
transfer function result. The reading error multiplied by one order of magnitude below the 
induced current error. 


c. Tides (Pawlowicz et al, 2002) 

Using the error output of t_tide and averaging the errors of the 
constituents, for one month, we have an uncertainty of around 0.020m/s. If we use 3 days 
the uncertainty is almost the same, 0.022m/s. 

d. Growing Surface Waves (Ardhuin et al, 2003) 

This phenomenon may “induce stresses opposed to the wind stress for 
wave growth stages that may represent up to 10% of the wind stress for short fetches.” 
For a typical current, and according our transfer function, this represents an uncertainty 
of about 0.02m/s. 


e. Eddy Surface Velocity 

In this study, a typical eddy speed variation over 24h is considered to be at 
least 0.02m/s. Eddies are highly variable but usually typical of certain locations. This 
uncertainty depends on the physical characteristics and if the eddy is influencing 
movement near the core or the out skirts. Observing the relative path is also important. 

f Others 

Other kinds of possible interferences were considered, but not included; 
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(1) Stokes Drift. For gravity waves: 

U ^ cok (15) 

where a is the amplitude, co the frequeney and k the wave number (Radko, 2007, Oeean 
Dynamies elass notes NFS). 

This yields for a Im amplitude, 10s period, 10m length wave an 
indueed veloeity of 0.39 m/s, wtiieli means an order of magnitude similar to the total 
RMS error. Studies show that it may be opposed and eaneelled by a Eulerian flow (Xu 
and Bowen, 1993; Monismith et ah, 2006). Works are still in progress to evaluate the 
total waves i-dueed flow. 

(2) Inertial Movements. These play a role but depend on the 
loeal flow oharaeteristies. They do not seem to play a major influenoe in the Monterey 
buoy area eireulation. Its surfaee eurrents may be over O.I m/s, again the order of our 
RMS error. We have to take into aeeount that if using yesterday’s eurrent for tomorrow, 
this error will be minimized. If we use tide treatment in a zone of important inertial 
motions, eare should be taken before averaging for the mean flow. 

(3) Other Non Linearities. (Eddies, larger periodieity waves, 
ete.) All of these influenees exist and are a part of the flow. The use of persistency 
already takes into account enough of its magnitude. 

The total sum of the uncertainties is 0.14m/s. Once these values 
were quantified in a conservative way, this means that the predictable accuracy is not 
expected to be better than the order of our actual results. 

2. Abrupt Changes 

The plotting of RMS error against significant wave height (SWH) and current 
modulus (following Garfield et ah, 2007) is almost intuitive and, by visual inspection, 
points out possible error causes. Wave, wind and current variations are shown as time 
series in Eigure 35, together with the RMS forecast current error time series. Whenever 
an abrupt change exists in the current speed, RMS seems to go up. The wind acts against 
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this rise, decreasing the RMS each time it couples the change in the current. Wave height 
seems not to have a significant consistent impact. This impression is clearly supported in 
the zoomed time series presented in Figure 36. 


RMS forecast with persistence and wind 



Figure 35. Full RMS plot versus SWH, current and wind. (Monterey station, Jul 06 to 
Dec 07.) 


60 

































































RMS forecast with tides treatment and wind 



RMS forecast with persistence and wind 



Figure 36. Zoom of RMS versus SWH, eurrent and wind - High error seems associated 
with large current variation. Top with tides treatment. Bottom with 
persistence. (Monterey station.) 
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IV. CONCLUSIONS 


This study analyzed long records of HF radar-derived surface currents and 
mooring-based wind observations for the region offshore Monterey Bay, California. The 
growth of HF radar instrumentation and the resulting surface current mapping data has 
enabled a number of scientific and public safety capabilities. In many instances, 
particularly with respect to search and rescue or spill mitigation applications, surface 
current forecasts are critical in addition to real time observations. To make accurate 
forecasts requires knowledge of background flow patterns plus variable flow patterns 
related to tides and wind-driven currents. The focus here has been on created and testing 
models for wind-driven surface currents that are specific to the near-surface (~0.5m) 
depths of the HF radar observations. 

The Monterey Bay region contains a wealth of surface observations from both HF 
radar and meteorological buoys. One site, in particular, surrounding the NOAA/NDBC 
buoy 46042 location has provided an eighteen-month, continuous set of currents and 
winds. Those data were used to develop an empirical transfer function between wind and 
HF radar-derived current. Shorter data sets from other buoy sites both inside and outside 
Monterey Bay were used to confirm results from the buoy 46042 location. In addition to 
the forward problem of creating a best-fit transfer function from wind to current, the 
recommended parameters were adjusted slightly by determining which values provided 
the most accurate 24-hour forecast of HF radar-derived surface current at the buoy 46042 
location. The results of these wind versus surface current analyses from the 1.5 year 
dataset from five different stations lead to the following conclusions and requirements. 

A. THE WIND INFLUENCE 

The wind induced current is on the order of 2% of the total wind speed, which is 
consistent with historical values. Despite this magnitude, the influence of introducing an 
explicit wind-driven current correction to the existing mean plus tidal current forecast 
models leads to around 10% of increased accuracy in the forecasts. The best results are 
obtained with a 50“ clockwise rotation, and with a magnitude = 0.01 St/j+0.018. 
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Although this study has shown a clear benefit to ineluding a wind-driven correetion, this 
eorreetion is still below the level of the foreeast uneertainties that derives from other 
uneorreeted faetors. 

B. RESULT CAPABILITIES 

The overall result, ineluding the wind-driven eurrent eorreetion, shows that HF 
radar-derived surfaee eurrent mapping data are able to foreeast surfaee eurrents in the 
next 24h with errors of about 0.15 +/- 0.1 m/s with 95% eonfidenee. When signifieant 
ehanges in eurrent happen without the wind influenee, approximately one day of high 
RMS foreeast error is expeeted. The largest and most persistent error events for the test 
period from July 06 to Deeember 07 led to a four-day RMS error over 0.2 m/s. 

C. NEXT STEPS 

The goal here has been to develop improvements to the short-term surfaee eurrent 
foreeast model, with primary foeus on ineluded the effeets of predieted wind ehanges. 
Several other faetors that are expeeted to influenee the aeeuraey of 24-hour surfaee 
eurrent foreeasts were pointed out in Chapter III. The results here show that it will be 
neeessary to understand and ineorporate more of those infiuenees if improvements 
beyond the 10% level aehieved here are desired. The qualitative assessment of error 
sourees presented in the last ehapter, eombined with previous work on HF radar surfaee 
eurrent validation, suggests that errors in the measurements themselves may still be the 
biggest limitation to improved foreeasts. 

In addition to work that is ongoing to improve HF radar surfaee eurrents and error 

traeking, this study points to the need to begin to ineorporate the most important aspeet of 

HF radar observations into the foreeast proeess. That is, spatial information from the 

surfaee eurrent mapping data should be ineorporated. At this time, all information about 

the sub tidal, tidal, and wind-driven surfaee eurrents that goes into the foreeast models 

derives from time series information from a single loeation. HF radar, on the other hand, 

provides a map of surfaee eurrents every hour. Observations show those maps eontain a 

good deal of horizontal variation in the surfaee eurrents, whieh is the fundamental reason 

that a network of HF radar systems is neeessary for the transport applieations supported 
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here compared with, for example, a single, moored current meter. The next generation 
forecast models should account for advective effects. The use of spatial integration 
between one calculated point and its immediate neighbors may allow gross errors to 
emerge more clearly. 
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